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We have employed a mesh experiment for back-illuminated (BI) CCDs. BI CCDs 
possess the same structure to those of FI CCDs. Since X-ray photons enter from the back 
surface of the CCD, a primary charge cloud is formed far from the electrodes. The primary 
charge cloud expands through diffusion process until it reaches the potential well that is 
just below the electrodes. Therefore, the diffusion time for the charge cloud produced is 
longer than that in the FI CCD, resulting a larger charge cloud shape expected. 

The mesh experiment enables us to specify the X-ray point of interaction with a subpixel 
resolution. We then have measured a charge cloud shape produced in the BI CCD. We 
found that there are two components of the charge cloud shape having different size: 
a narrow component and a broad component. The size of the narrow component is 
2.8 — 5.7 /im in unit of a standard deviation and strongly depends on the attenuation 
length in Si of incident X-rays. The shorter the attenuation length of X-rays is, the larger 
the charge cloud becomes. This result is qualitatively consistent with a diffusion model 
inside the CCD. On the other hand, the size of the broad component is roughly constant 
of ~ 13 fim and does not depend on X-ray energies. Judging from the design value of 
the CCD and the fraction of each component, we conclude that the narrow component is 
originated in the depletion region whereas the broad component is in the field-free region. 

KEYWORDS: charge-coupled device, X-ray event, split event, subpixel spatial resolution, 
mesh experiment 
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1. Introduction 

A charge-coupled device (CCD) is widely used in an optical imaging as well as an X-ray 
imaging. When an X-ray photon is photoabsorbed in the CCD, it generates a primary 
charge cloud consisting of thousands of electrons. The size of the primary charge cloud is 
73x (E/2.3 keV) 1 ' 75 nm where E is energy of X-ray photon in unit of keV. 1 ) The primary 
charge cloud expands through diffusion process until they reach the potential well of the 
CCD pixel. The final charge cloud after diffusion process is collected into several pixels 
forming various types of event pattern depending on how they split. 

We have developed a new technique "mesh experiment" which enables us to restrict 
the X-ray point of interaction with a subpixel resolution. 2 ) Hiraga et al. 3 ^ investigated 
the event pattern produced and measured the charge cloud shape produced in the front- 
illuminated (FI) CCD for the first time. They found that a charge cloud shape could 
be well represented by a Gaussian function. They also obtained the standard deviation, 
a, of the final charge cloud to be 0.7 ~ 1.5 ^m for 1.5 ~ 4.5 keV X-rays by using a FI 
CCD. Based on this experiment, they confirmed that there are three parameters tightly 
coupled together; 4 ) an X-ray point of interaction within a pixel, an X-ray event pattern 
and a charge cloud shape. Any two parameters can determine the third one. The event 
pattern is quite easily noticed while the charge cloud shape is difficult to be measured. 
Currently, they can separately measure it only by using the mesh experiment. Therefore, 
they can determine the X-ray point of interaction with much better spatial resolution 
than the pixel size. They obtained the position resolution of 0.7/im using the pixel size of 
12 fim if an X-ray photon became a split pixel event. However, due to its relatively small 
size of the charge cloud, the fraction of split pixel events whose point of interaction can 
be improved is less than 10 % of the total events. 

In this paper, we applied the mesh experiment to the back-illuminated (BI) CCD. 
Figure 1 shows the cross-section and the potential profile of FI CCDs. In the case of 
FI CCDs, since X-ray photons enter from the front side where there are electrodes. A 
major part of X-rays are photo-absorbed close to electrodes, resulting a final charge cloud 
to be relatively small. Figure 2 shows the cross-section and the potential profile of BI 
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CCDs. In the case of BI CCDs, a substrate of CCD, shown in Fig 1 (a), is etched from 
the back toward the depletion layer (thinning process) in order to achieve high quantum 
efficiency for blue lights and low energy X-rays. After the thinning process, almost all 
of the substrate has been removed. The thickness of the BI CCD employed is ~ 20 fim. 
Since X-ray photons enter from the back surface of the CCD shown in Fig 2, a primary 
charge cloud is formed far from the electrodes. The primary charge cloud expands through 
diffusion process until they reach the potential well which is just below the electrodes. 
Therefore, the diffusion time for the charge cloud generated by X-rays is longer than those 
in the FI CCD, resulting that a larger charge cloud shape is expected. 

Low energy X-rays, however, are photo-absorbed in the field-free region of BI CCDs. 
The primary charge cloud expands due to a diffusion process toward the electrodes as 
well as the back surface. The charge traveled to the back surface will be lost due to a 
recombination process. In the production of BI CCDs, accumulation is the most significant 
process that gives an internal potential to the back surface of CCDs and repels the signal 
charge to the electrodes. 5 ^ In the case of the BI CCD employed, the accumulation process 
is performed by ion implantation. 6 ) As shown in Fig 2 (b), the potential at close to the 
back surface is locally lower than that at the inner region caused by the accumulation 
process which enables us to collect the charge generated close the back surface. 

2. Experiment and Results 

The basic idea and the algorithm of data reconstruction for the mesh experiment are 
given in Tsunemi et al. 2 ^ We performed mesh experiment employing a BI CCD (S7170) 
fabricated by Hamamatsu Photonics K.K. The BI CCD has a two-phase electrode and 
full-frame transfer type. It possesses 512x512 active pixels with a size of 24 /im square. 
The mesh made of gold has small holes of 2.1 /im diameter with spaced 48 /im apart. The 
pitch of mesh holes is twice larger than that of the CCD pixel. The mesh is placed just 
5 mm above the CCD. We employed the 21m long X-ray beam line in our laboratory. 
We used the X-ray generator of Ultra-X18, fabricated by RIGAKU, with Al, Mo, and Ti 
targets in order to irradiate the characteristic emission line of each target: Al-K (1.5 keV), 
Mo-L (2.3 keV), and Ti-K (4.5 keV). We applied the voltage of 6, 5, and 7kV for Al, Mo, 
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and Ti target, respectively. We used a mechanical shutter made of stainless steel having 
thickness of 60 /im to control the X-ray beam so that the pile-up would not become a 
serious effect. We drove the BI CCD with our newly-developed system named E-NA 
system. 7 ^ The CCD analog data were processed by an integration-type circuit 8 -' and 
digital data were acquired with a VME system. We controlled the operating temperature 
of the CCD to be — 100 °C during the experiment. 

Figure 3 shows X-ray spectra of single-pixel events obtained with three targets em- 
ployed. The readout noise of our system is ~20e~rms. Each spectrum shows the char- 
acteristic emission line superimposed on the continuum emission. We should note that 
the peak channel of Mo-L emission line is close to that of Al-K emission line whereas the 
energy is quite different. 

We selected events between two dashed lines shown in Fig 3 in the subsequent analysis. 

2. 1 Event pattern 

Figure 4 shows a part of a raw image obtained with single-pixel events. We can see a 
clear moire pattern. Based on the moire pattern obtained, we can determine the geomet- 
rical relation between the mesh and the CCD. Details of the method are given in Tsunemi 
et a/. 9 ) 

Once we find their relations, we can determine the X-ray point of interaction within a 
pixel for all X-ray events. We thus rearranged each X-ray event with a subpixel resolution 
and investigated the distribution of single-pixel event, two-pixel split event, three- and 
four-pixel split event and more extended events as a function of the point of interaction 
inside a pixel. Figure 5 shows their X-ray intensity maps in 3x3 representative pixels 
(RP) for Mo-L X-rays. Each black dot corresponds to an individual X-ray event. When 
X-ray photons enter the center of a pixel, a major part of them produce single-pixel 
events shown in (a). When the X-ray point of interaction becomes close to the horizontal 
or vertical boundary of each pixel, the horizontally or vertically two-pixel split events are 
formed shown in (b) or (c). When an X-ray photon enters the corner of each pixel, it 
produces three- or four-pixel split event shown in (d). This tendency is consistent with 
that obtained for the FI CCD. 10 ) However, we found that the regions forming the three- 
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and four-pixel split events become larger and extend into inner region than FI results (e.g. 
Figure 7 of Tsunemi et al. 9 ^) whereas the pixel size of the BI CCD employed is 24 fim 
that is two times larger than that of the FI CCD. This suggests that a charge cloud size 
produced in the BI CCD is much larger than that of the FI CCD. This fact is clarified in 
Fig 5 (e) since a charge cloud spreads into more than four pixels even if an X-ray enters 
the center of a pixel. 

2.2 Absorption structure of BI CCDs 

As shown in Fig 1, X-ray photons enter from the electrodes in the case of FI CCDs. The 
thickness of electrodes is not uniform and some part of electrodes is overlapped, resulting 
the detection efficiency within a pixel to be far from uniform. 2 ' n ~ 13 ) It makes it difficult 
to make the response matrix of FI CCDs. 

As to the BI CCDs, on the other hand, X-ray photons enter from the back surface 
where the thickness of the field-free region is much thinner than those of FI CCDs shown in 
Fig 2. The potential created by the accumulation process also drives the charge generated 
in the field-free region toward the electrodes, resulting higher detection efficiency for low 
energy X-rays. Since the thickness of the back surface of CCDs is uniform, we can expect 
a uniform detection efficiencies within pixels for BI CCDs. Figure 6 shows the X-ray 
intensity map extracted from all X-ray events in 3x3 RPs for Mo-L X-rays. There is no 
significant variation in each image, demonstrating the uniform detection efficiency of the 
BI CCD. We also confirmed uniformity for Al-K and Ti-K X-rays. The non-uniformities 
in detection efficiency are ~ 2.5, 2.0, and 7.2% in a standard deviation for Al-K, Mo-L, 
and Ti-K X-rays, respectively. This uniformity is a large advantageous point of BI CCDs 
over FI CCDs. We should note that the uniformity of Ti-K X-rays is lower than those of 
other X-rays. It must be due to the fact that the mean attenuation length of Ti-K X-rays 
is long enough for X-ray point of interaction to be affected by electrodes. 
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2.3 Amount of charge collected in a center pixel 

When X-rays enter the CCD at the position of (X p , Y p ), the output of the n th pixel of 
the CCD, P n (X p ,Y p ), is written as 

P n (X p , Y p ) = dX dY f(X — X p , Y — Y p ) 

where (X n , X n+ i, Y n , Y n+ i) denotes the boundary of the n t h pixel and f(X, Y) is the 
charge cloud shape just before the collection by the potential well. The actual output of 
n th pixel, D n (X in , Y in ), is a convolution between P n (X p , Y p ) and the mesh hole, which can 
be described by 



D n (X in , Y in ) — J dXp J dYp P n (X p , Y p ) H(X in — X p , Y in — Y p ) (2.1) 

= dX dY / <g> H(X - X m , Y - Y m ) (2.2) 

where (X in , Y in ) is the center of mesh hole, H(X, Y) is a typical hole shape of the mesh, 
and / (g> H represents the convolution between / and H. 

We can experimentally measure D n (X in , Y in ) by using the RP generated with all X-ray 
events. Figure 7 show the surface profiles and projected profiles of D n in 3x3 pixel regions 
where the center pixel is n t h pixel for Mo-L, Al-K, and Ti-K X-rays. 

Hiraga et al. 3 ^ measured D n for the FI CCD for the first time and found a flat-top 
feature around the center of n th pixel. This is due to the fact that the charge collected 
in the n t h pixel is almost constant, forming single-pixel events, when a mesh hole is well 
within the pixel. D n obtained with Ti-K for the BI CCD has relatively flat top at the 
center of pixel. On the other hand, D n shows a sharp structure even near the center of 
pixel in Mo-L and Al-K X-rays, suggesting that a charge splits into neighboring pixels 
even if X-rays enter near the center of the pixel. We can again expect that a charge cloud 
size produced in the BI CCD is much larger than that of the FI CCD. It is interesting 
to note that a flat-top feature obtained only with Ti-K X-rays suggests that a size of the 
charge cloud of Ti-K X-rays is smaller than those obtained with Mo-L and Al-K, which 
is in contrast to those obtained with FI CCDs. 10 ) 
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2.4 Charge cloud shape 

Hiraga et al. 3 ^ calculated a charge cloud shape by differentiating D n (we hereafter 
refereed to the differential method). As shown in Fig 7, however, we cannot obtain a 
charge cloud shape with the differential method since D n does not show a flat-top feature 
for the BI CCD. We then employed a different approach to obtain a charge cloud shape. 

As written in eq 2.1, D n is an integration of convolution between a charge cloud shape 
and a mesh hole within a CCD pixel. We directly calculate their convolution to reproduce 
D n for each energy assuming a charge cloud shape to be a two-dimensional axial-symmetric 
Gaussian function 3 ^ (the integral method, hereafter). We should note that we have applied 
the integral method to the FI data and obtained consistent results to those obtained by 
the differential method. We found, however, a difficulty to reproduce D n with the integral 
method. The problem is that our method cannot reproduce a sharp component and a tail 
component, simultaneously. Such a tail component is detected for all energies whereas it 
was not detected with FI CCD. 3 ) The presence of an tail component suggests a large size 
of a charge cloud which must be much larger than that forming a sharp component. As 
shown in Fig 5 (e), the events spreading more than 4 pixel is uniformly distributed within 
a pixel. These events should account for such tail components. 

We thus introduce an extra component to fit the data. We employed two components 
in the charge cloud shape having different size in both axes: a narrow component for 
the sharp structure and a broad component for the tail structure. In this way, we found 
that the two-component model in the integral method well reproduced D n obtained. The 
model functions of D n are shown in Fig 8. Table I summarizes the charge cloud sizes 
of two components for all energies. Figure 9 shows the size of the charge cloud for both 
components as a function of the X-ray attenuation length in Si. 

The broad component of the charge cloud is 2.5 — 5 times larger than that of the nar- 
row component. The size of the narrow component decreases as the attenuation length 
increases whereas that of the broad component is relatively constant of ~ 13 /im. Com- 
paring the FI results, 10 ) the sizes for the narrow and the broad component are 2—8 times 
larger and an order of magnitude larger, respectively. The number of X-ray events of the 
broad component is ~ 3 times larger than that of the narrow component. 
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3. Discussion 

3.1 Non-linearity of pulse height versus energy 

As shown in Fig 3, the pulse height of each characteristic X-ray line does not linearly 
increase with energy. Figure 10 shows their relationship between the pulse height and the 
incident X-ray energy. The pulse height of Mo-L is ~ 30 % lower than that expected from 
the linear relationship based on Al-K and Ti-K X-rays as shown by a solid line. This is 
quite a different characteristics of the BI CCD to the FI CCD. 

The attenuation length of Mo-L X-rays is only 2.2 /im that is the shortest among X-rays 
employed. Therefore, most of Mo-L X-rays are absorbed close to the accumulation region, 
where the density of impurities is an order of magnitude larger than that in the depletion 
region, leading to the recombination between impurities and signal charge. We thus 
suppose that the non-linearity of pulse height versus energy is caused by an incomplete 
collection of the signal charge produced in the accumulation or some part of the field-free 
region. 

3.2 Origin of the two components of charge cloud 

We found that there were two components in a charge cloud shape generated by the 
BI CCD. A narrow components has a er of 2.8 — 5.7/im which strongly depends on the 
incident X-ray energy. The size of a broad component is relatively constant of ~ 13 /im 
having less dependence on the X-ray energy whereas such a component is not appeared 
in the FI CCD. 3 ) 

As to BI CCDs, X-rays entering from the back surface of the CCD form a primary 
charge cloud far from the electrodes. The primary charge cloud expands through diffusion 
process until they reach the potential well which is just below the electrodes. Therefore, 
a diffusion time for the charge cloud produced by low energy X-rays is longer than that 
by high energy X-rays, resulting a larger size of charge cloud shape for low energy X-rays. 
This relation is consistent with that obtained for a narrow component. However, we found 
that the size of the broad component is almost independent of the incident X-ray energy. 

As shown in Figs 1 (b) and 2 (b), the potential profiles of FI and BI CCDs are different 
in the field-free region. A charge generated in the field-free region can be collected only 
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through the diffusion process in FI CCDs. On the other hand, the potential created by the 
accumulation drives the charge generated in the field-free region toward electrodes in the 
BI CCD. Therefore, a charge generated both in the depletion region and in the field-free 
region can be collected in the BI CCD. With taking into account the difference in the 
charge cloud size for both components, we suppose that an X-ray photo-absorbed in the 
depletion region can be detected as a narrow component whereas an X-ray photo-absorbed 
in the field- free region can be detected as a broad component. 

Miyata et al. 8 ^ measured the thickness of the depletion depth of FI CCD having similar 
wafer as that of the BI CCD employed to be ~ 4//m. Since the thickness of the BI 
CCD employed is 20 jum, thickness of the field-free region is ~ 16 /im. This value is in 
good agreement with the charge cloud size of the broad component if the charge cloud 
expands spherically in the field-free region. Moreover, as shown in Table I, the fraction 
of the broad component obtained is ~3 times larger than that of the narrow component, 
which is also consistent with the ratio of thickness between the field-free region and the 
depletion region. These facts strongly support that a broad component is originated in 
the field-free region. 

3.3 Asymmetry of the charge cloud shape 

The charge cloud size obtained with the BI CCD is shown in Table I. Both for the 
narrow and broad components, the charge cloud shapes are almost point symmetric. The 
difference in both directions is less than 5%. This is also noticed that D n obtained for 
three energies shown in Fig 7 show a quite point-symmetric shape. 

Hiraga et a/. 3 ) pointed out that a charge cloud shape obtained for the FI CCD shows 
an asymmetric shape. The size in unit of the standard deviation for the x-direction 
was 1.4 /im whereas that for the y-direction was 0.7 — 0.8 fim. Tsunemi et al. 14 ^ further 
investigated the possible explanations about the asymmetry of the charge cloud shape. 
They conclude that the asymmetry is due not to the experimental setup, such as the 
asymmetric mesh hole, but to the asymmetry of the electric field inside the CCD. Similar 
asymmetry has obtained for some other FI CCDs. 10,11 ) Therefore, the asymmetric shape 
of the charge cloud is commonly found in FI CCDs. 
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On the other hand, we obtained a symmetric shape of the charge cloud generated inside 
the BI CCD for all three X-rays. The back surface of CCDs are uniform and is normally 
applied to be zero voltage. At the front side, there are four electrodes per pixel in our 
case and they are slightly overlapped to each other as shown in Fig 1 (a). Moreover, 
at the front side, there is another gate structure, "channel stop", which forms the pixel 
boundary along the y-direction. All these gate structures generate the electric field not to 
be uniform near the front electrodes as pointed out by Tsunemi et al. 1 ^ The symmetry of 
the charge cloud obtained with BI CCD suggests that a charge cloud must not be affected 
by the electrodes significantly because it is generated far from the electrodes. 

4. Conclusion 

We have performed the mesh experiment with the BI CCD and measured the charge 
cloud size for the first time. We found that there are two components of the charge cloud 
shape to reproduce the data obtained all for Mo-L, Al-K and Ti X-rays: a narrow com- 
ponent and a broad component. The narrow component possesses a a of 2.8 — 5.7 /xm 
which strongly depends on the attenuation lengths of incident X-ray energies. The broad 
component, on the other hand, possesses much larger a than that of the narrow compo- 
nent. It is ~ 13 (Km and almost independent of the attenuation length of X-rays. The 
accumulation technique to fabricate BI CCDs is not employed in the FI CCDs and the 
charge produced even in the field-free region can be collected toward electrodes. We thus 
suppose that the narrow and the broad components are originated in the depletion region 
and the field-free region, respectively. This hypothesis is consistent with both by the 
design value of the CCD and by the fraction of X-ray events for both components. 

We have found that the linearity between the incident X-ray energies and the pulse 
height is not good. The pulse height of Mo-L X-rays is ~ 30 % lower than that expected. 
This might be explained by the fact that the effect of the recombination with the impurities 
is significant in the accumulation region since the attenuation length of Mo-L is shortest 
among those X-rays employed. 
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Figure captions 

Fig. I. (a) Cross-section and (b) potential profile of a FI CCD. 
Fig. 2. Same as Fig 1 but for a BI CCD. 

Fig. 3. X-ray spectra obtained with targets of (a) Al, (b) Mo, and (c) Ti. Each charac- 
teristic emission lines are identified. The dotted lines show the ranges employed for 
analysis. 

Fig. 4. A part of the raw image (200 x 200 pixels) extracted from the single-pixel events 

of Mo-L X-rays. Each dot shows an X-ray event. A clear Moire pattern can be seen. 

Since the pitch of the mesh hole is two times that of the CCD pixel, we find X-ray 

events appeared in every 2x2 pixels. 
Fig. 5. X-ray intensity map of (a) single-pixel events, (b) horizontally two-pixel split 

events, (c) vertically two-pixel split events in 3x3 RPs, (d) three- and four-pixel split 

events, and (e) more extended events with the subpixel resolution. The dotted lines 

show the pixel boundary. 
Fig. 6. X-ray intensity map extracted from all X-ray events in 3x3 RPs for Mo-L. The 

uniform detection efficiency can be found for BI CCD. Projected profiles onto X- and 

Y-axes are also shown. The dotted lines show the pixel boundary. 
Fig. 7. Amount of charge collected in the n t h pixel, D n , is shown in 3x3 pixels for (a) 

Mo-L, (b) Al-K, and (c) Ti-K X-rays. Upper panel of each figure shows the 3-d picture 

of D n and projected profiles are shown in lower panel. 
Fig. 8. Model of D n for (a) Mo-L, (b) Al-K, and (c) Ti-K X-rays. 
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Fig. 9. Charge cloud size as a function of the X-ray attenuation length in Si for (a) the 

narrow component and (b) the broad component. 
Fig. 10. Relationship between the pulse height and incident X-ray energy. The solid line 

shows the linear relationship between them to reproduce Ti-K X-rays. 
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Table I. Charge cloud shape obtained with the BI CCD for Mo-L, Al-K, and Ti-K X-rays. 





Mo-L 


Al-K 


Ti-K 


Energy [keV] 


2.3 


1.5 


4.5 


Attenuation length [/zm] 


2.2 


7.9 


13.5 


Narrow component 


&nx 


5.3 ±0.3 


4.3 ±0.2 


2.91 ±0.08 


&ny 


5.8 ±0.3 


4.4 ±0.2 


2.99 ±0.09 


Snx 


5.3 ±0.3 


4.2 ±0.1 


2.83 ±0.08 




5.7 ±0.3 


4.3 ±0.2 


2.91 ±0.09 


Broad component 




13.4 ±0.3 


13.2 ±0.3 


11.9 ±0.1 


&wy 


13.9 ±0.3 


13.8 ±0.3 


12.5 ±0.1 


£>wx 


13.4 ±0.3 


13.2 ±0.3 


11.9 ±0.1 


Swy 


13.9 ±0.3 


13.8 ±0.3 


12.5 ±0.1 


Fraction of narrow component 




0.25 ±0.02 


0.32 ±0.01 


0.32 ±0.01 



15 



This figure "fl.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/physics/0206021vl 



This figure "f2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/physics/0206021vl 



This figure "f3a.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/physics/0206021vl 



This figure "f3b.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/physics/0206021vl 



This figure "f3c.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/physics/0206021vl 



This figure "f4.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/physics/0206021vl 



This figure "f5a.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/physics/0206021vl 



This figure "f5b.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/physics/0206021vl 



This figure "f5c.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/physics/0206021vl 



This figure "f5d.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/physics/0206021vl 



This figure "f5e.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/physics/0206021vl 



This figure "f6.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/physics/0206021vl 



This figure "f7.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/physics/0206021vl 



This figure "f8.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/physics/0206021vl 



This figure "f9a.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/physics/0206021vl 



This figure "f9b.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/physics/0206021vl 



This figure "flO.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/physics/0206021vl 



